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Pion form factor using domain wall valence and asqtad sea quarks 
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We compute the pion electromagnetic form factor in a hybrid calculation with domain wall valence quarks and 
improved staggered (asqtad) sea quarks. This method can easily be extended to p — + 77r transition form factors. 
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The pion electromagnetic form factor is con- 
sidered a good observable for studying the onset 
with increasing energy of the perturbative QCD 
regime for exclusive processes. As the pion is the 
lightest and simplest hadron, a perturbative de- 
scription is believed to be valid at lower energy 
scales than predictions for heavier and more com- 
plicated hadrons like the nueleon p. 

A pseudoscalar particle has only a single elec- 
tromagnetic form factor, F{Q^), where is the 
four-momentum transfer, and at = 0, this 
form factor is the electric charge of the parti- 
cle, F{Q^ = 0) = 1; the magnetic form factor 
vanishes. The experimentally observed behavior 
of the form factor at small momentum transfer 
is well described by the vector meson dominance 
(VMD) hypothesis pTllj 
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The current experimental situation is presented 
in Fig.[n ,5.6.7.8.,9) . 

What is surprising is that VMD with only the 
lightest resonance (mp = mvMo) appears to de- 
scribe all the existing data, even up to scales of 
> IGeV^. In contrast, at very high momen- 
tum transfer, we expect the data to approach the 
perturbative behavior [10111112] 
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Higher order perturbative calculations of the hard 
contribution to the form factor [13114115] do not 
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Figure 1. Summary of experimental data for the 
pion electromagnetic form factor. Shaded regions 
are expected sensitivities of future experiments. 

vary significantly from this value, as seen in 
Fig. n At the largest energy scale where reliable 
experimental measurements have so far been ob- 
tained, around ~ 2 GeV^, the data are 100% 
larger than this pQCD asymptotic prediction. 

Early lattice calculations validated the vec- 
tor meson dominance hypothesis at low [161 
[TTj . Recent lattice results [18I19I2()I21| . includ- 
ing some of our own preliminary results [22I23| . 
have somewhat extended the range of momentum 
transfer, up to 2 GeV^, and the results remain 
consistent with VMD and the experimental data. 

For our exploration of the pion form factor with 
unquenched gauge configurations, we performed 
a hybrid calculation using MILC Nf = 2+1 
and Nf = 3 configurations in 20'^x64 volumes. 
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Table 1 

Simulation details for domain wall fcrmion calcu- 
lations on 20^x64 dynamical MILC asqtad lat- 



tices at a 



1.6 GcV. 
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Figure 2. Pion (left) and rho meson (right) dis- 
persion relation vs. continuum (upper) and lattice 
(lower) expectations curves. 



generated with staggered asqtad sea auarksPl]. 
and domain wall valence quarks with domain wall 
height TTiQ = 1.7 and extent Lg ~ 16 of the extra 
dimension The MILC configurations were 

HYP blocked before valence propagators were 
computed, otherwise the residual chiral symme- 
try breaking would have been unacccptably large. 
Dirichlet boundary conditions were imposed 32 
timeslices apart. Thus, opposite halves of the 
lattices were used on alternate configurations in 
the Monte Carlo sequence in order to reduce au- 
tocorrelations. A detailed study of the physical 
properties of light hadrons composed of staggered 
quarks computed on these lattices has recently 
been completed |27|. Details of the observables we 
computed are available in our earlier references 
22 23j. 

To compute the pion form factor at large mo- 
mentum transfer, we must first understand the 
pion dispersion relation when the momenta are 
relatively large. In the continuum limit, we ex- 
pect proper relativistic behaviour 



El{v) = v' + El{0). 



(3) 



From the study of free lattice bosons, we can de- 
fine lattice equivalents of the continuum momen- 
tum and energy 



E = 2 sinh {E/2) = 2 sin (p,/2) . 



(4) 



Another possibility is that the pion dispersion re- 
lation will follow that of a free lattice boson 



Em 
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Eqs. (O and ((SJ) differ significantly only when the 
lattice momenta are large. 

In Fig. 12 we have plotted against each disper- 
sion relation and we see that the data clearly favor 
the lattice dispersion relation at large momenta. 
This gives us confidence that we can use Eq. © 
to determine £'7r(p) at high momenta from fits to 
low momentum correlators which have much less 
statistical noise. 

Using the lattice dispersion relation and the ra- 
tio method described in our previous work, we 
have computed in Fig. 13 the pion form factor 
for two dynamical pion masses. Fits of the data 
to the monopole form of Eq. are also shown, 
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Figure 3. Pion electromagnetic form factor for 
fixed sink momentum = (1, 0, 0) computed by 
the ratio method |22I23| and imposing the lattice 
dispersion relation, Eq. Shaded regions are 
jackknife error bands for VMD fit. 

where shaded regions correspond to jackknife er- 
ror bands and the central values are given in the 
legend. The data in Fig. |21 were computed with 
pseudoscalar pion sink operator fixed at momen- 
tum = (1,0,0). We were unable to compute 
statistically significant form factors with the same 
pseudoscalar sink operator for p^ = (1,1,0) by 
the ratio method, apparently due to the poor 
overlap of the sink operator with the higher mo- 
mentum state. We expect that an axial vector 
sink operator would be a better choice for future 
calculations. A comparison of the ratio method 
with the fitting method, also described in j22l23| . 
should appear shortly |2H1- 
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